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Superconductor strip with transport current: Magneto-optical study of current 

distribution and its relaxation 

A. V. Bobyl 1 ' 2 , D. V. Shantsev 1 * 2 , Y. M. Galperin 1 ' 2 , T. H. Johansen 1 *!, M. Baziljevich 1 , and S. F. Karmancnko 3 
1 Department of Physics, University of Oslo, P. 0. Box 1048 Blindern, 0316 Oslo, Norway 
2 A. F. Ioffe Physico- Technical Institute, Polytekhnicheskaya 26, St.Petersburg 194021, Russia 
3 Electrotechnical University, Prof. Popov Str. 5, St. Petersburg, 197376, Russia 

The dynamics of magnetic flux distributions across a YBa2Cu307_5 strip carrying transport 
current is measured using magneto-optical imaging at 20 K. The current is applied in pulses of 
40-5000 ms duration and magnitude close to the critical one, 5.5 A. During the pulse some extra 
flux usually penetrates the strip, so the local field increases in magnitude. When the strip is 
initially penetrated by flux, the local field either increases or decreases depending both on the 
spatial coordinate and the current magnitude. Meanwhile, the current density always tends to 
redistribute more uniformly. Despite the relaxation, all distributions remain qualitatively similar to 
the Bean model predictions. 



I. INTRODUCTION 

In many technological applications like high-T c 
bolometers and electrical power cables the superconduc- 
tors are subjected to pulses of transport current. If the 
current is small, the magnetic flux behavior in super- 
conductors with strong bulk pinning is described by the 
critical-state model (CSM).p According to the model, the 
magnetic flux is frozen, and the current density j is ev- 
erywhere equal to (or less than) its critical value, j c . 
Flux and current distributions in superconductors carry- 
ing a transport current were measured in several previous 
works. Among them are. results obtained using-naagneto- 
optical (MO) imaging, Bi3 Hall micro-probes ,E2Hl3 and 
detection based on THz radiationO 

If the transport current / exceeds its critical value I c 
the CSM description becomes invalid. In this case the 
flux starts to move through the superconductor resulting 
in a non-zero electric field across the sample, i. e., in a 
finite resistance. In reality, the transition from frozen to 
moving flux lattice is rather smooth, since for any value 
of I some regions of the superconductor carry the cur- 
rent density j ~ j c . For / < I c these regions occupy 
a substantial fraction of the sample, and the magnetic 
flux distribution becomes highly non-stationary. This is 
known to result in a number of dynamic effects such as 
voltage dependence xui the current sweep rate, and resis- 
tance relaxation.EJ~EZl Direct spatially-resolved measure- 
ments of the non-stationary flux distribution could be 
very useful for analysis of these effects. Such measure- 
ments have already proved quite efficient for the applied 
field case, namely, to access the spectral distribution of 
the activation energies in YBa2Cu307_a filmsJlS and to 
detect transinet vortex phases in Bi-2212 crystals.Ej In 
the present paper we use the MO imaging technique to 
study the dynamics of flux and current distributions in- 
duced by a pulsed transport current / up to I c . The 
experiments are performed on a thin strip of YBa2Cu3- 
07_5 and cover the time interval 40-5000 ms after onset 
of the pulse. 

It turns out that the most convenient regime to study 



the dynamic effects is when the strip is subjected to both 
a transport current and an externally applied field. Then 
one can create states where the CSM predictscj a re- 
gion with j — j c covering more than a half of the strip 
width, and therefore the deviations from the CSM are 
more easily detected. Such -a- situation has so far been 
studied only in a few workspElllj and without a system- 
atic comparison to the CSM and with no attention to 
dynamic effects. Therefore, we have chosen the following 
two magnetic states for the present study: (i) a flux-free 
strip subjected to a transport current - the reference case, 
and (ii) a strip in the remanent state after a large field 
pulse plus a transport current. 

The paper is organized as follows. The experimental 
method including the field-to-current inversion scheme is 
described in Sec. ||. Section III presents comparison of 
the current density distributions for different I to the 
CSM predictions. The relaxation of the flux and current 
distributions at fixed I is analysed in Sec. IV. 



II. EXPERIMENTAL PROCEDURE 

Films of YBa2Cu3p,7_5 (YBCO) were grown by dc 
magnetron sputteringciJ on LaA103 substrate. X-ray and 
Raman spectroscopy analysis confirmed that the films 
were c-axis oriented and of a high structural perfection. 
A microbridge with dimensions 500 x 100 x 0.2 /im 3 was 
formed from the film by a standard lithography proce- 
dure. Several Au wires of 50 [mi diameter were attached 
to each Au/Ag contact pad using thermocompression. 

The imaging system is based on the Faraday rotation 
of polarized light illuminating an MO-active indicator 
film that we mount directly on top of the superconduc- 
tor's surface. The rotated Faraday angle varies locally 
with the value of the perpendicular magnetic field, and 
through an analyzer in an optical microscope one can di- 
rectly visualize and quantify the field distribution across 
the covered sample area. As a Faraday-active indica- 
tor we use a Biidoped ferrite garnet film with in-plane 
magnetizationEa The indicator film was deposited to a 
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thickness of 5 /jm by liquid phase epitaxy on a gadolin- 
ium gallium garnet substrate. A thin layer of aluminum 
is evaporated onto the film allowing incident light to be 
reflected, thus providing double Faraday rotation of the 
light beam. The images were recorded with an 8-bit Ko- 
dak DCS 420 CCD camera and transferred to a computer 
for processing. The MO imaging at low temperatures 
was performed by mounting the superconductor/MO- 
indicator on the cold finger of a continuous He-flow cryo- 
stat with an optical window (Microstat, Oxford). 

The flux distribution in a sample biased with large 
transport current is subjected to fast relaxation. There- 
fore, it was important to minimize the camera exposure 
time. On the other hand, very short exposures lead to 
poor-quality MO images, and as a compromise, 33 ms 
exposure time was chosen. The sample was biased with 
current pulses of 40 ms duration, which were synchro- 
nized with the camera recording. One image was taken 
during the pulse, and another image after the current 
was switched off. Experiments were also performed to 
study the relaxation of flux distributions during current 
pulses. For this purpose we used pulse durations up to 
5000 ms. The MO images were then recorded at different 
times during the pulse. 

As the bridge thickness is much less than its width, the 
theoretical results for the thin strip geometrjeja can be 
used. In our experiments, the bridge thickness is also 
of the order of the London penetration depth, hence its 
magnetic properties are fully characterized by the two- 
dimensional flux distribution at the surface. 
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FIG. 1. Superconducting strip with transport current 

The current density distribution in the strip cross sec- 
tion was calculated from the measured flux density distri- 
bution using the ID inversion scheme proposed in Ref. |24|. 
Consider a long thin superconducting strip with edges lo- 
cated at x = ±w, the y-axis pointing along the strip, and 
the z-axis normal to the strip plane, see Fig. ^. The trans- 
port current flows in y-direction as does the screening 
current due to the magnetic field, B a , applied along the 
z-axis. The perpendicular magnetic field at the height h 
above a superconductor strip, i. e. at — w < x < w, is 
related to the current density distribution by the Biot- 
Savart's law 



Here J(x) = J j(x, z) dz is the sheet current, j(x,z) is 
the current density, and the integral is calculated over 
the strip thickness Z, where Z <C w. Inversion of Eq. (|l|) 
yieldso 
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This expression is written for discrete coordinates x = 
nA with integer n, the ratio h/A is denoted as d. The 
current distributions J(x) calculated from this expres- 
sion are already reasonably good, see, e. g., Refs. pip. 
However, since the high-fc components are suppressed by 
the inversion procedure, the accuracy of the calculated 
J{x) turns out to be poor near the sample edge where 
the current drops abruptly to zero. 

In the present work we develop an improved inversion 
scheme using an iteration procedure. Suppose we know 
Jn{x) calculated on the Nth iteration step. We then set 
Jn(x) equal to zero outside the superconductor, and sub- 
stitute the new current distribution into Eq. (Q) to get the 
flux distribution B^(x). The latter will be slightly differ- 
ent from the measured distribution B(x). The difference 
B(x) — Bn(x) is then substituted into Eq. (||) to calcu- 
late 5 Jn (x). We then set Jn+i(x) = Jn(x)+SJn(x) and 
start the next iteration. After each iteration step Bn(x) 
approaches closer to B{x), and the current near the edge 
drops more abruptly. It is usually sufficient to perform 
5-10 iterations. The final current distribution exhibits an 
abrupt jump at the sample edges as can be seen from the 
figures below and a test example discussed in Appendix 
A. In the calculations we also took into account effect of 
the component B x on the MO indicator sensitivityO 



III. PULSES WITH DIFFERENT MAGNITUDE 

A. Transport current 

The MO images for an initially flux-free strip carrying 
a transport current were taken at 20 K. The image for 
the highest current, / = 4.9 A, is shown in Fig. |2[ The 
image was taken with slightly uncrossed polarizer and 
analyzer. This provides that the flux density of different 
polarity shows up as different gray level. Consequently, 
the MO images can be directly converted into a flux den- 
sity distribution using a careful calibration curve. The 
flux distributions were calculated for the area marked by 
the black rectangle in Fig. [2] and then averaged along the 
strip. 
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Six profiles of the flux density B{x) across the strip 
for different currents are shown in Fig. |^(a). The flux 
density has opposite signs at the left and the right parts 
of the strip, and \B\ is maximal near the edges. The 
corresponding current density profiles, Fig. ||(b), are cal- 
culated by the inversion scheme, using h = 12/im. The 
experimental current profiles reproduce well the main 
features of the Bean-model result, Fig. ||(c), where the 
current density is always minimal at the center of the 
strip. On the other hand, the model predicts constant 
J(x) = J c near the edges. Instead, experimentally we 
find maxima in J(x) their magnitude and position be- 
ing dependent on /. Tbis-behavior is in agreement with 
flux creep simulations,c§lJ hence these deviations from 
the CSM can be related to dynamic effects. 




FIG. 2. MO image of YBCO strip carrying transport cur- 
rent 4.9 A. The image is obtained with uncrossed polarizers 
so that the positive and negative flux density correspond to 
bright and dark regions, respectively. The black rectangle 
shows the image area used for calculation of the current den- 
sity distribution. 

When the transport current I is switched off, the mag- 
netic flux in the inner part of the strip remains trapped. 
The return field of this trapped flux re-magnetizes the 
edge regions of the strip and flux of opposite sign pen- 
etrates an outer rim. Flux density profiles for such 
a remanent state taken across the strip are shown in 
Fig. |](a) . The central part of these profiles is similar to 
that of the corresponding profiles in the current-carrying 
state, Fig. ||(a). However, typical values of \B\ are much 
smaller, in agreement with the Bean model, therefore the 
scale of y-axis is different. 

The current density profiles are shown in Fig. [|(b) 
along with the Bean- model profiles, panel (c). As the 
net current must be zero in any cross section, current 



lines should form closed loops. In the center, the cur- 
rent keeps flowing in the same direction as the transport 
current, while near the edges it now flows backwards. 
This result is nicely reproduced in our experiment. More- 
over, in agreement with the Bean model we find that (i) 
the current density is maximal near the edges, (ii) the 
boundary between positive and negative J shifts inwards 
as / increases. Nevertheless, like in the current-carrying 
state, we don't find an agreement with the theory near 
the edges where J(x) should be always equal to J c . 
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FIG. 3. Initially flux-free YBCO strip carrying transport 
current 7=1.8, 2.8, 3.5, 4.3, 4.6, and 4.9 A. (a,b) experimental 
flux and current density profiles determined from MO images 
like in Fig. ^, (c) the CSM current density profiles, Eq. (|l|), 
with I c = 5.5 A. The flux density data were obtained from 
MO images of the strip, and the plot B(x) is the average 
profile in a 300 /jm wide band across the strip. 

The flux dynamics is expected to proceed much faster 
in the regions where the current density is close to the 



3 



critical one. It is exactly those edgy regions where the 
deviations from the Bean behavior are found. However, 
in the above experiment, the size of these regions was 
relatively small, 0.6 — 0.9w < \x\ < w. Thus, the results 
could possibly be affected by inaccuracy of the inversion 
method or small defects near the strip edge. For this 
reason we performed another experiment where the re- 
gion of j ~ j c is much more extended, and unambiguous 
manifestations of the dynamic effects are found. 
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FIG. 5. YBCO strip, initially in fully-penetrated state, car- 
rying transport current 1=0, 3.5, 4.3, 4.6, and 4.9 A. (a,b) ex- 
perimental flux and current density profiles, (c) the CSM cur- 
rent density profiles, Eq. Q), with I c = 5.5 A. Inset in (b) is 
a blow-up of the dashed rectangle. 



FIG. 4. Remanent state of initially flux-free YBCO strip 
after a current pulse with 7=2.8, 3.5, 4.3, 4.6, and 4.9 A. 
(a,b) experimental flux and current density profiles, (c) the 
CSM current density profiles, Eq. (□) with I c = 5.5 A. 



B. Remanent Field + Transport current 

In the second experiment, a large magnetic field was 
applied to the strip and subsequently removed. As a re- 
sult, some flux remain trapped in the strip, and the flux 
profile has a conventional triangular shape with maxi- 
mum at the center, see the "I = 0" curve in Fig. ||(a). 
Small regions near the edges, |x| > are filled with 

return field of the opposite polarity. The strip was then 
biased with transport current I, and flux density dis- 
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tributions were measured in the current-carrying state. 
B(x) profiles across the strip for different / are shown 
in Fig. ||(a). As the current increases, the central peak 
diminishes, and almost disappears at the maximal cur- 
rent (solid line). Meanwhile, another peak at the left edge 
shows up and grows with /. At intermediate currents one 
can clearly see two peaks in B(x) profile. This feature, 
was predicted long ago within the Bean model analysis ,t3 
but has not so far been observed in experiment. 




FIG. 6. Remanent state of initially fully-penetrated YBCO 
strip after a current pulse 1=0, 1.8, 3.5, 4.3, 4.6, and 4.9 A. 
(a,b) experimental flux and current density profiles, (c) the 
CSM current density profiles, Eq. (^), with I c = 5.5 A. 



The current density profiles are presented in Fig. ||(b) 
along with the Bean model profiles, panel (c). According 
to the model, the current distribution in the left half is ac- 
tually a copy of the distribution for the "pure" transport 
current state, Fig. ||(c), though squeezed in x-direction. 



Thus, J(x) profile has a deep minimum at the middle of 
the left half, which is well reproduced in the experimen- 
tal profiles. However, we find a clear disagreement with 
the theory in the right half of the strip. There the Bean 
model predicts a flat J{x) unaffected by the transport 
current. In experiment, J(x), though being quite flat ap- 
pears dependent on /. This is emphasized in the inset, 
where one can see a monotonous decrease of J(x) level 
as the current increases. Such a behavior contradicts not 
only the Bean model but also common intuition. This in- 
teresting observation is discussed below in the framework 
of relaxation effects. 

The remanent flux and current distributions after the 
current pulse are shown in Fig. ^. One can clearly see 
a current-induced suppression and shift-pf the peak in 
B(x), aa predicted by the Bean model fB and observed 
earlier .113 However, the current density profile is too com- 
plicated to make any solid conclusions. We mention only 
two main features of the Bean profiles found also exper- 
imentally: (i) a peak in J(x) at small negative x, (ii) a 
concave shape of the profile in the left, and convex in 
right half of the strip. 



IV. RELAXATION DURING THE PULSE 

In order to understand the reason for the observed de- 
viations of the measured B- and J-profiles from the CSM 
predictions we analyze the effect of pulse duration. It was 
found that the flux distribution does not remain fixed 
during the current pulse, even though the current magni- 
tude is kept constant. The flux is subjected to relaxation 
which is especially fast for larger currents. 

The relaxation effects are known to be noticeable in 
high-temperature superconductors even in the absence 
of transport current due to large flux creep rates .E3 The 
direction in which the relaxation proceeds can be always 
well understood. In a conventional experiment, the flux 
trapped in a superconductor kept at zero applied field 
monotonously decreases with time. This means that the 
flux density decreases in amplitude throughout the sam- 
ple. If a sample is kept in constant non-zero applied field 
H ai then the flux distribution relaxes towards a uniform 
H(r) — H a . Remarkably for a thin sample in a perpen- 
dicular field there coexist regions with H > H a , where H 
decreases dufi-to, creep, and regions with H < H a , where 
H increases. E3E3 They are separated by a so-called neu- 
tral line, where the flux density remains constant. 

In the presence of transport current the relaxation of 
flux and current density profiles is even more compli- 
cated as illustrated by Figs. 0, and ||. Flux distribu- 
tions were measured at different times during the pulse 
starting from t = 40 ms up to 5000 ms (t = corre- 
sponds to switching current on) . Fig. ^ shows the results 
for the case of initially flux-free strip with I = 4.9 A, 
like in Fig. ||. Only the first (t — 40 ms) and the last 
(t = 1600 ms) B(x) profiles are presented. The two cor- 
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responding J(x) profiles are indicated by thick solid and 
thick dashed line respectively at the bottom panel. 
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FIG. 7. Relaxation of flux and current density distribution 
in an initially flux-free YBCO strip during the current pulse, 
/ = 4.9 A. The profiles correspond to different times: 40, 100, 
200, 300, 450, 800, 1200, and 1600 ms; arrows show the time 
direction. The absolute flux density increases throughout the 
sample. The current tends to redistribute more uniformly: 
some extra current moves from the edges to the central re- 
gion. 



We note that, in contrast to the applied field case, 
the relaxation leads to increasing the field amplitude 
throughout the superconductor. This means that more 
and more flux enters the sample as the time passes. In 
terms of J-profile it means that some extra current re- 
distributes from the edges to the center. It can be also 
interpreted as a tendency for J-distribution to become 
more uniform. 



Shown in Fig. g are results for the combined "remanent 
field + transport current" state, with / = 3.5 and 4.9 A. 
Here the relaxation of B{x) does not follow a simple rule. 
The flux density may increase or decrease depending on 
the position in the sample, and on the current ampli- 
tude. Moreover, a closer look shows that the direction 
of relaxation may change as the relaxation proceeds, like 
in Fig. ||(a) ati» —0.8w. However, the changes in J- 
profile have the same tendency as for the flux-free strip 
subjected to transport current: the current tends to dis- 
tribute more uniformly. 
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FIG. 8. Relaxation of flux and current density distribution 
during the current pulse for the "remanent field + current" 
state. Top: / = 4.9 A, times: 40, 100, 200, 800 and 1600 ms. 
Bottom: I = 3.5 A, times: 40, 400, 800, 1600, and 5000 ms. 
Arrows show the time direction. The current tends to redis- 
tribute more uniformly: some extra current moves from the 
high-current right half to the low-current left half. 



In Refs. 29 3(], the relaxation of the local flux den- 
sity measured by a Hall sensor was described using an 
"effective" CSM with a time- dependent critical current 
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density. This description is basically in agreement with 
our results. Indeed, when J c decays with time, the CSM 
predicts more and more uniform J(x)-distribution un- 
til at some point one comes to a completely uniform 
J(x) = J c {t). The rate of the J c decay can be estimated 
from the plateau in J(cc)-profiles at x > in Fig. [| We 
find that from t = 40 to 5000 ms the J c reduction was 
» 10%. 

We have also performed computer simulations of flux 
creep based on the Maxwell equation, and Bio-Savart's 
law for a thin strip, Eq. (|l]), ancLiagarithmic current 
dependence of the pinning energyEo The direction in 
which the relaxation of flux and current proceeds was al- 
ways in agreement with our experimental results. The 
same tendency has been also found in simulations re- 
ported in Ref. [T^ . 

However, the creep simulations do not explain the de- 
crease of J with transport current increase illustrated in 
the inset of Fig. B. One can see from comparison of the 
top and bottom panels of Fig. || that this apparent de- 
crease is related to a faster relaxation of J-distribution for 
larger /. Consequently, the relaxation rate of J c turns out 
to be /-dependent. The explanation of this effect remains 
an open question. One can expect that at larger currents 
a slow permanent flow of vortices across the strip sets 
in, probably along a few channels. This permanent flow 
does not change the global flux distribution and cannot 
be resolved magneto-optically. However, it keeps vortices 
in motion which should reduce J c and assist relaxation 
processes. 



V. CONCLUSION 

We have studied the dynamics of flux and current den- 
sity distributions in a YBCO strip subjected to a pulse of 
transport current at 20 K. Relaxation of the distributions 
during the pulse is much different from the conventional 
relaxation of trapped flux in the absence of current. The 
local flux density may decrease or increase depending on 
the position in the sample, current amplitude, and mag- 
netic history. However, the relaxation of J-distribution 
always has the same tendency: the current tends to re- 
distribute more uniformly. Despite the relaxation, most 
flux and current distributions remain in a qualitative 
agreement with the Bean model predictions and can be 
described by a Bean model with a time-dependent J c . 
However the behavior for the current-carrying strip in 
initially fully-penetrated state is not fully understood. 
This behavior is suggested to be associated with a flux- 
flow induced reduction of J c . 
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APPENDIX A: THE INVERSION METHOD 



We demonstrate here advantages of the proposed it- 
eration inversion method for calculation of the current 
density profile J{x) across the strip from the measured 
field profile. For this purpose we use a model example 
of a strip with Bean- model J{x) given by Eq. (|l|), with 
a = 0.5w. The corresponding flux density profile B(x) at 
the height h above the strip is given by Eq. (|l|) . To model 
our experimental situation we assume that B(x) is mea- 
sured in the range —Aw < x < Aw at fixed x separated 
by A = 0.02w, and that h = 0.2w. 

First, we restore the current density profile using the 
standard procedure based on Eq. (Q). Straightforward 
calculations where d — h/A = 10 produce highly oscil- 
lating J(x) due to contribution of higher harmonics, see 
pi| . The best results are obtained by the choice d = 2, 
which means that only B-values at x separated by h/2 
are used. This J{x) profile shown by the dashed line in 
Fig. H reproduces qualitatively the initial J-profile (dot- 
ted line). There are, however, large distortions near the 
edges, and an error in the overall current magnitude. 

The solid line shows the J-profile obtained by the iter- 
ation procedure described in Sec. II. C, for 5 iterations. It 
is clear from the figure that the iteration method is more 
successful in restoring the correct J-profile. It should be 
noted that the exact form of Eq. @ is not essential for 
the iteration procedure. It can be modified if it will help 
improve convergence of B^(x) towards B{x). The best 
convergence is found when parameter d is equal to unity 
or less. 
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FIG. 9. Current density distributions obtained by itera- 
tive inversion procedure (solid line), and standard procedure 
(dashed line) for a Bean-model field distribution. The dotted 
line is Bean-model current distribution. 



7 



The advantage of the iteration procedure is explained 
by the fact that the inversion formula, Eq. (§), is used 
to calculate only S Jn (x) instead of the whole J(x) . Nor- 
mally, 5Jn(x) is much smaller and smoother than Jn(x). 
Thus, an inevitable error intrinsic to any inversion pro- 
cedure is minimized. The only assumption used by this 
procedure is that current cannot flow outside the su- 
perconductor sample, which is physically well-grounded. 
The location of the sample edges can be determined with 



high accuracy from a simple optical image. 

APPENDIX B: BEAN-MODEL RESULTS 

Here we list some Bean-model expressions for current 
density distributions in a thin superconducting strip ob- 
tained in Refs. and ^3|. For an initially flux- free strip 
— w < x < w, Fig. [I], carrying transport current I 



M. = I 1 arctan ■ N < « , (i) 

Jc 1, a < \x\ < w 

where a = Wy/l — (I/I c ) 2 , and I c = 2w J c is the critical current. This distribution is shown in Fig. ||(c). In the 
remanent state, after the transport current is switched off, 



J(x) _ j | [arctan (J$=g) - 2 arctan (J&g) 



-a < x < a, 

(2) 

• ? . 1 1 - | arctan (^/|£f£) , a < \x\ < b 



where a — Wy/l — (I/I c ) 2 , b — wy/l — (I/2I C ) 2 . At b < \x\ < w the current density is equal to — J c , see Fig. [|(c). 

When a large magnetic field was applied to the strip and subsequently removed, the current density J{x) = — J c 
and J{x) — J c in the left and the right half of the strip, respectively. If the strip is subsequently biased with transport 
current /, then 



J ( x ) _ [ ~ arctan (^ giggE ^) - 1, ~w/2 - a < x < -w/2 + a, 

Jc I 1, —w < x < —w/2 — a or — w/2 + a<x<w 

see Fig. |5](c). In the remanent state after the current pulse, see Fig. ^|(c), 

a < x < 6, 



(3) 



M = / "I " f -tan (^fe£) + | arctan a t{Z^ ) 



2 



Jc I — 1 b < \x\ < w 



(4) 
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